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^ ; ABSTRACT 

m 

The cumulene carbenes are important components of hydrocarbon chemistry 
in low mass star forming cores. Here we report the first astronomical detection 
of the long chain cumulene carbene H2C6 in the interstellar cloud TMC1, from 
observations of two of its rotational transitions: Jr,k' — 7i,7 ~~ ► 61^ at 18.8 GHz 
and 8^8 — > 7ij at 21.5 GHz, using NASA's Deep Space Network 70 m antenna 
at Goldstone, California. In addition we also observed the shorter cumulene 
carbene, H2C4 at the same position. The fractional abundance of H 2 Cg relative 
to H 2 is about 4.7 x 10 -11 and H 2 C 4 is about 1.1 x 10~ 9 . The abundance of H 2 C 6 
is in fairly good agreement with gas phase chemical models for young molecular 
cloud cores, but the abundance of H 2 C4 is significantly larger than predicted. 
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1. INTRODUCTION 



The detection of new cumulene carbenes is important for determining the hydrocarbon 
chemistry in low mass star forming cores (cf. Bettens & Herbst 1996). The cumulene 
carbene chains are also of interest in astronomy for their possible role as carriers of the 
diffuse interstellar bands (McCarthy et al. 1997). Here we report the first astronomical 
detection of H 2 Cg (hexapentaenylidene) in the dark cloud TMC1. The identification of this 
highly reactive isomer of triacetylene (HCeH) was made possible by the recent laboratory 
detection of the rotational spectrum of both H2C5 and H2C6 (McCarthy et al. 1997), and 
the high sensitivity of NASA's Deep Space Network antennas at cm wavelengths. 

The shorter cumulene carbenes propadienylidene H2C3 (Cernicharo et al. 1991, 
Kawaguchi et al. 1991) and butatrienylidene H2C4 (Kawaguchi et al.) were also observed 
at the same position as H2C6, allowing comparison with abundances predicted by chemical 
models. Both H2C4 and H 2 C6 have linear carbon backbones and C2*, symmetry and, 
because of the off-axis equivalent H- atoms, they both have ortho (K = ±1, odd) and para 
(K = 0, even) rotational ladders between which both radiative and collisional transitions 
are forbidden. In cold molecular clouds we expect 3/4 of the population to be in the lowest 
lying ortho (K = ±1) ladders, and emission from these to be slightly stronger than from 
the para (K = 0) ladder. 

In addition to H2C6, several other complex molecules were observed in TMC1, for the 
purpose of clarifying hydrocarbon chemistry in these star forming cores. We find that the 
carbon chain molecules vary in abundance among the three principal velocity components. 
The distribution of the cumulene carbenes is similar to those of other complex carbon 
compounds, such as the cyanopolyynes HC2 n +iN and the carbon chain radicals C„H, but 
differ from the cyclopropenylidene ring, C-C3H2. The observed ratios of abundances in 
TMC1 are compared with gas phase chemical models. 



2. OBSERVATIONS 

The observations were made with NASA's Deep Space Network (DSN) 70 meter 
antenna at Goldstone, CA during October and November of 1996. In addition to H2C6, 
five other molecules (HC9N, C5H, H2C3, H2C4, and 0C3H2) were observed as part of an 
ongoing effort to study the carbon chemistry in TMC1. The relevant transitions and line 
frequencies are listed in Table 1. 

The observations were made with a broad band (17.5 - 26 GHz) cooled HEMT 
receiver with typical system temperatures of 50 - 65 K and the 2 million channel Wide 
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Band Spectrometer Analyzer (WBSA). The spectra were smoothed to resolutions in the 
range of 0.03 to 0.15 km s _1 with a total spectral range of 20 to 40 MHz (cf. Langer et 
al. 1995). All observations, except that of the 18.4 GHz transition of H 2 C6, were made 
using position switching. In this mode the spectra were obtained by taking the difference 
between successive ON and OFF source integrations. To suppress baseline irregularities of 
instrumental and atmospheric origin, the right ascension of the OFF source position was 
also chosen to be 6 minutes in time to the west or east of the ON source position such 
that each ON-OFF pair had the same mean elevation. The presence of several stronger 
lines in the passband of the 21.5 GHz line of H 2 C@ and the 20.8 GHz line of H2C3, allowed 
us to check our calibration during each observing session. Following the detection of the 
21.5 GHz line of H 2 Cg, a frequency switching option was implemented for the DSN-WBSA 
system. We used this mode for the observations of the 18.8 GHz transition of H 2 C6 using 
an LO offset of 0.5 MHz between the signal and reference frequencies. 

The detections reported here were made toward TMC1 at a position [a(1950) = 
4 h 38 m 41 s and 5(1950) = 25°35'40"] corresponding to the strongest HC 7 N and HC 9 N 
emission observed at t>LSR — 5.8 km s _1 in our DSN 70 m high spectral resolution maps 
(Velusamy et al. 1997). The ON source integration times for the H 2 C6 lines were 23 and 18 
hours for the 21.5 and 18.8 GHz transitions, respectively, and yielded a corresponding rms 
of 2.5 niK and 3.5 mK at spectral resolutions of 0.14 and 0.12 km s _1 . 



3. RESULTS 

In Figure la we show the spectra for the two transitions of H 2 C6, and those of C5H 
and HC9N. The H 2 C6 21.5 GHz transition, HC9N, and C5H were observed simultaneously 
in the same passband. The spectrum of the 2o, 2 — > lo,i transition of H2C4 at 17.9 GHz in 
TMC1 is shown in Figure lb, along with the lines of H2C3, CqH, and 0C3H2. The H2C4 
and H 2 C6 spectra peak at ~ 6 km s _1 with antenna temperatures T^ = 180 mK and 9 
mK, respectively. Previous high spectral resolution studies of carbon chain molecules in 
TMC1, such as CCS and HC7N (Langer et al. 1995), revealed that there are three velocity 
components at 5.7, 5.9, and 6.1 km s _1 with very narrow linewidths ~ 0.15 to 0.20 km s _1 . 
Our detection of H 2 C@ in TMC1 supports the suggestion of a tentative assignment of a 
weak unblended and a blended 3 mm line in IRC+10216 (Guelin et al. 1997). 

In Figure 1 it appears that the strongest H2C4 and H 2 C6 emission is at 5.9 km s _1 . 
The slight asymmetry in the line profiles is most likely due to a second weaker (partially 
blended) component at 6.1 km s _1 , exactly what is seen in the velocity structure of C5H 
and probably C@H. There is no evidence in H 2 C4 and H 2 Cg for a component at 5.7 km s _1 . 
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In contrast, C-C3H2 and its linear isomer H2C3 show strong emission at 5.7 and 5.9 km s -1 , 
and only slightly less emission at 6.1 km s _1 . The presence of strong 0C3H2 and H2C3 
emission at 5.7 km s" 1 and the apparent absence of H 2 C6 or H2C4 emission at that velocity 
indicates that there is a significant chemical inhomogeneity in the hydrocarbon chemistry 
in TMC1. Furthermore, the 0C3H2 emission is much stronger than that of H2C3 (Ta = 2.0 
K versus 0.12 K). 



4. DISCUSSION 

To estimate the abundances of H 2 C6 and H2C4 in TMC1 we adopted a simple model in 
which we treat each rotational ladder as a separate linear rotor. The excitation of the ortho 
and para states can be treated separately, because the radiative and collisional transitions 
between them are negligible. For purposes of estimating abundances, the excitation of 
the ortho ladder can be treated as two equivalent linear rotors. The energy levels of the 
K = ±1 ladders are nearly equal and can be estimated within 1% by treating them as a 
linear rotor with Bq = (B + C)/2 (see McCarthy et at). Collisional rate coefficients for 
these species have not been calculated; we approximate their J-dependent values by using 
those derived for the linear chain molecule HC3N (Green and Chapman 1978). We assume 
that the rotational population of H 2 C6 has a normal ortho:para ratio of 3:1. Furthermore, 
reactive collisions with protonating ions can produce H2C6H 4 " which, upon decomposition 
to H2C6 + H, will preserve the nuclear spins (Herbst, 1977). 

H 2 C6 is highly polar, with a calculated dipole moment of 6.2 D (Maluendes & McLean 
1992). The J = 8 ortho levels lie ~ 5 K above the ortho ground state (J = 1) and these 
levels have an Einstein A coefficient of about 2 x 10~ 6 s~ x . This yields a critical density 
n(li2) ~5x 10 3 cm -3 , implying that the lower rotational levels are easily thermalized at 
the density in the cores of TMC1. Figure 2 shows an excitation calculation for the K — — 1 
ortho ladder of H 2 C6 and the K = para ladder of H2C4 (using a dipole moment of 4.5 D) 
for a range of fractional abundance gradients (in km s _1 pc _1 ). For convenience, calculations 
were done with a large velocity gradient (LVG) code which is essentially equivalent to 
an LTE calculation for optically thin lines, as appropriate for the weak emission seen in 
H2C4 and H 2 C6. We assumed typical conditions found in the TMC1 clumps: Tj, in = 10 K 
and n(H 2 ) = 10 4 cm~ 3 (Langer et al.). The velocity gradient is estimated from the line 
width, ~ 0.3 km s _1 , divided by the size of the emission region (~ 0.06 pc) estimated from 
our HC7N map. The total fractional abundances of H2C4 and H2C6, after accounting for 
molecules in the other ladders, a velocity gradient of 5 km s _1 pc -1 , and a beam efficiency of 
0.7, are about 1.1 x 10~ 9 and 4.7 x 10 -11 , respectively. The column density for H2C4 derived 
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by Kawaguchi et al. (1991) from several transitions at a nearby position corresponds to a 
fractional abundance of 8 x 10~ 10 (Bettens et al. 1995), which is in very good agreement 
with our value. We estimate that H2C4 is ~ 25 more abundant than H 2 C6. 

Recent gas phase chemical models have included reaction networks for calculating 
abundances of long carbon chains (Millar et al. 1997, Bettens, Lee & Herbst 1995) and very 
large hydrocarbons (Thaddeus 1994, and Bettens & Herbst 1996). Furthermore, laboratory 
data and theoretical calculations have begun to clarify the mechanisms for forming the 
simplest ring structure cyclopropynylidyne, 0C3H, via neutral carbon insertion into 
acetylene, C + C2H2 (Kaiser et al. 1996). To compare qualitatively our observations with 
chemical models we use the results of the gas phase model of Millar et al. (1997), which is 
based on one of the standard chemical data bases used for modeling cloud chemistry. It 
includes hydrocarbon production up to 8 or 9 carbons, but does not include the carbon 
insertion reactions. 

In Figure 3 we plot these model abundances for the carbon chains C2 ra H, C2 n +iH, 
H 2 C 2 n, H 2 C 2n +i, and HC 2n +iN, as well as a point for c-C 3 H 2 , for the physical conditions 
relevant to this core of TMC1. Figure 3 uses the early time results of Millar et al. as 
only this evolutionary stage has enough neutral carbon to produce a successful complex 
carbon chemistry with large hydrocarbon abundances. The late times, or steady state, 
solutions have too little carbon and complex carbon species to explain the observations. 
This result is not surprising as other tracers, such as CCS, confirm the ages of these dense 
core components to be < a few xlO 5 years (Kuiper et al. 1996, Velusamy et al. 1997). 
Note that the C n H radicals are predicted to be more abundant than the corresponding 
H 2 C n cumulene carbenes, as observed. However, the Millar et al. model also predicts that 
H 2 C6 is more abundant than H2C4 which is exactly the opposite from what is observed. 
The measured H 2 C6 abundance agrees within a factor of 2 with their model calculations, 
while the predicted H2C4 fractional abundance is too low by a factor ~ 25. On the other 
hand, the observed H2C4/H2C6 ratio is consistent with the slope of the abundance versus 
chain length in the C n H radicals. However, some caution is in order regarding absolute 
abundances, because this model does not distinguish between different isomers, such as 
triacetylene (HCeH) and H2C6. 

Bettens et al. explored the effects of various assumptions about the neutral-neutral 
reactions on the production of complex carbon molecules in interstellar clouds. Bettens 
& Herbst (1996) recommend using two of their chemical models which best explain 
the polyatomic species: the new standard model (NSM) and model 4 (M4). NSM is a 
modified version of the standard ion-molecule chemical scheme; in M4, rapid neutral-neutral 
reactions play a critical role. A key feature of M4 is that the reactions O + C n — ► 
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C„_i + CO are assumed negligible for n > 2. Model M4 of Bettens et al. provides the 
best agreement with our H 2 C4 observations, but the calculated abundance is still too 
high by a factor of five. The models appear to be sensitive to assumptions about the 
neutral-neutral reaction rate coefficients at low temperature. H2C6 and H2C4 are important 
diagnostics to discriminate among models of hydrocarbon chemistry in dense cores, however, 
more laboratory measurements and astronomical observations are needed to resolve the 
hydrocarbon chemistry in dense cores. 

We thank the staff of the DSN at Goldstone for their assistance in making these 
observations and Eric Herbst for useful comments. The research of the JPL group was 
conducted at the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. WDL would also like to 
thank the Smithsonian Astrophysical Observatory and the Center for Astrophysics for their 
hospitality during two visits while this work was in progress. 
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TABLE 1 



MOLECULAR TRANSITIONS OBSERVED IN TMC1 



Species 




Transition 


Frequency 
MHz 


Reference 


H2C3 


lo,i 


^0 ,o 


20792.640 


Lovas et al. 


c-C 3 H 2 


li,o 


->lo,i 


18343.143 


Pickett et al. 


H2C4 


2o,2 


-lo,i 


17863.810 


Tr avers et al. 


C 5 H 


J = 


9/2 -»• 7/2, F = 5 - 46 


21484.710 


Lovas 


C 5 H 


J = 


7/2 -»• 7/2, F = 4 - 36 


21485.262 


Lovas 


C 6 H 


J = 


15/2 — > 13/2,^ = 8-7/ 


20792.872 


Lovas 


H2C 6 


7l,7 


->• 61,6 


18802.235 


McCarthy et al. 


H2C 6 


81,8 


-7i, 7 


21488.255 


McCarthy et al. 


HC 9 N 


37- 


+ 36 


21498.181 


Pickett et al. 
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Fig. 1 — (a) DSN 70 m detections of H 2 C 6 at 21.5 and 18.8 GHz in TMC1. The spectra for 
H 2 C6 at 21.5 GHz, HCgN, and C5H, were taken within the same passband, while the 18.8 
GHz line of H 2 C6 was observed separately. The source position is a(1950) = 4 A 38 m 41 s and 
5(1950) = 25°35'40", with v LSR = 5.8 km s _1 . (b) The spectrum for H 2 C 4 towards TMC1 is 
shown along with those of CgH (2 hyperfine components), H 2 C3, and c-C3H 2 . The vertical 
line marks 5.9 km s _1 . 



Fig. 2. — The antenna temperatures for the ortho levels of H 2 Cq (solid lines) and para 
levels of H 2 C4 (dashed lines) calculated using a LVG excitation model and approximating 
the energy levels with a simple linear rotor (see text) as a function of J. The parameters are 
Tfci n = 10 K, n(H 2 ) = 10 4 cm -3 , and AV/AL = 1 km s" 1 pc _1 . Curves for three fractional 
abundances (in units of km s _1 pc _1 ) are shown for each molecule. The observed T^ for 
H 2 C4 (filled triangle) and H 2 Cg (filled circles) are marked in the figures along with their 
la error bars. Total fractional abundances need to include a factor for the ortho and para 
fractions, the velocity gradient, and antenna efficiency (see text). 



Fig. 3. — The fractional abundances of five carbon chains, C 2n H, C 2n+ iH, H 2 C 2n+ i, H 2 C 2n , 
and HC 2n+ iN as a function of number of carbons C n . These are taken from the early time 
solutions for the gas phase chemical model of Millar et al. for a cold cloud core with T^ in = 
10K and n(H 2 )=10 4 cm~ 3 . The observed fractional abundances of H 2 C4 and H 2 C6, assuming 
AV/AL = 5 km s _1 pc _1 , are marked by filled triangles. 



